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We present an overview of the temperature, field and angular dependence of the irreversible 
55 ' magnetization of non-magnetic borocarbides (Y;Lu)Ni2B2C. We show that nonlocal electrodynamics 

influences pinning via the unusual behavior of the shear modulus in non-hexagonal lattices. On top 
of that, we observe that the pinning force density F p exhibits a rich anisotropic behavior that sharply 
contrasts with its small mass anisotropy. When H_Lc, F p is much larger and has a quite different H 
dependence, indicating that other pinning mechanisms are present. 



I. INTRODUCTION 



Since early 1994, when the borocarbide family (RM2B2C, where R = rare earth) was discoveredETu a series of 
interesting results have been reported on these materials. The most remarkable feature is the coexistence of super- 
i conductivity and antiferromagnetism for R = Tm, Er, Ho and Dy. The high superconducting transition temperatures 
T c and the broad variation of the ratio Tn/T c (where Tm is theJNfeel temperature, ranging from 1.5 K to 10 K) make 
. this family particularly appropriate to explore that coexistenceou. 

Interesting properties are also found-ia the non-magnetic borocarbides, Rp^Y,Lu. These materials exhibit a non- 
30 and a non-linear field dependence!!^! 



^■—^ , exponential temperature dependenc and a non-linear field dependencet3 of the electronic pSoecific heat in the 
superconducting phase; an anomalous upper critical field i£;2 with an upward curvature near T c E3o and a foruufojd 
oscillation when the field is rotated within the basal planet^; and a square flux line lattice (FLL) at high-fieldst^rllll. 

Although these properties were initially taken as evidence of a non conventional pairing mechanismlla, both the 
unconventional H C 2 and the square FLL can also be satisfactorily explained by assuming a non-local electrodynamics. 
Traditionally, nonlocality was expected to be relevant only in materials with a Ginzburg-Landau parameter k = A/£ ~ 
1 and sufficiently clean to have an electronic mean free path t 3> £. However, the large vortex cores make theoretical 
analysis very difficult in that case. The borocarbides, in contrast, have intermediate k values (~ 10 — 20) that make 
O ■ core effects much smaller. This allows tha. non-local effects to be introduced perturbatively either in the London 
model, as dona-by Kogan, Gurevich et al.Ej, or in the Ginzburg-Landau theory, as shown by dc Wilde et al.Ej and 
Park and HuseEU. In both cases, the result is a coupling of the supercurrents tor-thc underlying crystal symmetry. 

The first success of the Kogan- Gurevich nonlocal model was the explanation^ of the deviations of the reversible 
magnetization M of Bi:2212 in the mixed state from the logarithmic dependence on magnetic field, M cx ln(H C 2/H), 
predicted by the London mode£3. Later on, Song et al.El found similar deviations in YM2B2C when H || c-axis, 
that could also be quantitatively accounted for by the model. Recently we have extended that study to all directions 
of H and showed that this generalization of the London Iheory provides a satisfactory complete description of the 
anisotropic M(H) with a self-consistent set of parameters^. 

In the local London model anisotropy is introduced via a second rank mass tensor my . In tetragonal materials such 
as borocasbides m a = nit, thus the properties should be isotropic in the basal plane. However, non-local corrections 
introduce^] a fourfold anisotropy as a function of the magnetic field orientation within the ah plane, refleptiug the 
square symmetry This 7r/2 periodicity was indeed observed in M in the mixed state of both YM2B2C0EJ and 



LuNi 2 B 2 CEa. 
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Kogan's model also predictedEa that two structural transitions in the FLL should occur in borocarbides for H | c, 
a first order reorientation transition between two rhombic lattices at a field Hi aud a second order transition from 
rhombic to square at H2 > Hi. Small angle neutron scattering (SANS) studiestLZl confirmed those predictions in 
YM2B2C. A jump in the apical angle (3 of the rhombic lattice, discontinuous within the resolution, occurs at Hi ~ 1 
to 1.25/cOe, aiuLthe lattice becomes square (/3 = 90°) at H2 ~ 1.25 to 1.5fcOe. According to a recent analysis by 
Knigavko et alHj, Hi really consists of two second order transitions taking place in a very narrow field range. 

Although the role of nonlocality on the equilibrium properties of the FLL is by now convincingly established, less 
is known about its effects on the nonequilibrium vortex response. Since nonmagnetic borocarbides exhibit a-very low 
critical current J c for H || c, pinning correlation volumes were expected to be large, as indeed observedEJ. Thus, 
the elastic properties of the FLL must play a key role in the pinning. As the shear modulus Cqq depends on pE\, 
and (3 undergoes a discontinuous jump at Hi, it is bound to happen that Cqq and therefore the pinning properties 
change abruptly at this fieldEH In other words, vortex pinning, which involves distortions from equilibrium vortex 
configurations, should he affected by the symmetry changes in the vortex lattice. 

Recently we showedE3 that, in YNI2B2C and for H || c, the reorientation transition at Hi induces a kink in the 
pinning force density F P (H). We also observed anisotropics in F p both between the c-axis and the basal plane (out- 
of-plane anisotropy) and within the plane (in-plane anisotropy). We found that F p for H_Lc is one order of magnitude 
larger than for H || c and has a quite different field dependence. We argued that this surprising behavior is unlikely to 
arise either from pinning by magnetic impurities or from non-local effects, although we recognized that the evidence 
supporting those claims was not conclusive. 

Here we report further studies of F p in non-magnetic borocarbides. We show that the kink in F p (H) is also 
visible in LuNi 2 B 2 C and coincides with the field Hi for this compound, thus confirming that it is a signature of 
this nonlocality-induced transition. We find that H\{T) slightly decreases as T increases, in contrast to H 2 (T). We 
also study the effect of Co-doping in Lu(Nii_2;Co x )2B2C. We observe that Hi decreases as the nonlocal effects are 
progressively suppressed by increasing x, in agreement with the T dependence. All the LuNi2_ ;c Co ;c B2C samples 
exhibit an enormous out-of-plane anisotropy. This unambiguously demonstrates that this anisotropy is due neither 
to the magnetic impurities (as those crystals have a density of impurities much smaller than the YNi2B2C crystal), 
nor to nonlocality (as it does not dissapear with increasing x). We also rule out the precence of surface barriers for 
H±c by performing minor hysterisis loops. 



II. EXPERIMENTAL DETAILS 



We report magnetization measurements, performed with a Quantum Design SQUID magnetometer, on YNi2B2C 
(Y-0), and Lu(Nii_ x Co x ) 2 B 2 C with x = (Lu-0), x = 1.5% (Lu-1.5) and x = 3% (Lu-3) single crystals. The 
dimensions, T c and estimated I ior_eacb .sample are summarized in Table I. The Y-0 crystal is the same one that was 
previously investigated in Refs.EjOc^OS. The normal state magnetization follows a Curie law that indicates the 
presence of a very dilute distribution of localized magnetic moments. It corresponds to a_rare-earth impurity content 
of 0.1 at. % relative to Y, probably due to contaminants in the Y starting materiao. The Lu(Nii_ a: Co a: )2B2C 
crystals, grown as described elsewherd±3, show a much weaker Curie tail at low_temperatures (T < lOOi-T) which 
might arise, for example, from a 0.001% magnetic impurities of Gd in the Lu siteEj. Isothermal magnetization loops 
in the superconducting mixed-state were measured, and the critical current density J c was then calculated using the 
Bean's critical state modelc-ffi3. 



III. RESULTS AND DISCUSIONS 
A. Reorientation phase transition 



Figure 0(a) shows the pinning force density F p — | J c x B| for the Y-0 sample as a function of the applied field for 



y with increasing H, but above a field 
we showed that this "kink" in F p (H) 



H || c-axis at several T. We observe that at low fields F p {H) decreases stror 
H* ~ 1.2fcOe the field dependence becomes much weaker. In a recent work 
at H * is a signature of the reorientation phase transition in the FLL. We based this claim in several facts. First, the 
position of the kink for H || c coincides with the value of Hi. Second, Jl* is rather insensitive to the field orientation, 
in agreement with the behavior ofrHi observed in SANS experimentsEH. Finally, Hi/H c i is an increasing function of 
T, as predicted by Knigavko et al.E3. 

TABLE I. volume, thickness, superconducting critical temperature and electronic mean free path for the investigated samples. 
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sample V(mm :i ) t(rnm) T C (K) 1(A) 

Y - 2l 05 153 300 

Lu- 2.5 0.3 15.7 270 

Lu- 1.5 1.2 0.4 14.9 100 

Lu-3 0.44 0.2 14.1 70 
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We now expand the previous-study and analyze the behavior of the field H* when the nonlocality is reduced. In the 
Kogan-Gurevich descriptiontJ'Ej the strength of the nonlocal perturbations is parametrized by a new characteristic 
distance, the nonlocality radius p(T, £): the weaker the nonlocality effects, the smaller p(T, £). This means that p(T, £) 
decreases with increasing T or decreasing t. Thus, if this model is correct, an increase in T or a decrease in £ should 
produce qualitatively similar effects on the nonlocality- induced transitions H\ and H2- 




H (kOe) 

FIG. 1. Field dependence of the pinning force density for H || c-axis for (a) Y-0 sample at several temperatures and (b) 
Lu-0, Lu-1.5 and Lu-3 samples at T = 5K. 

A typical way to reduce £ is to introduce impurities. JX has been shown that the increase of x in L^Nii-^Co^E^C 
crystals decreases £ without significantly increasing J c u. Recent resultSnConfirmed that the analogy between T and £ 
is indeed valid in this material in the case of i?2- First, Gammel et al.ES showed that H2 increases as £ is reduced by 
increasing x. Later on, Eskildsen et al.E3 showed that H2(T) also rises as T increases. 

If H* is indeed a signature of Hi, we should also find a correlation between its T and £ dependencies. In Figure |2| 
we plot H*(T), together with H C 2(T). We observe that H* is almost constant at low T and, unlike the rhombic to 
square second order transition H^iT), it slightly decreases with T at higher temperatures. 
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FIG. 2. Temperature dependence of the upper critical field and the structural reorientation transition field for the Y-0 sample 
when H || c. 

To determine the I dependence of H* we performed measurements on the Lu-0, Lu-1.5 and Lu-3 samples at T = 5K 
for H || c (estimates of £ are given in Table I). The results are shown in Figure |l](b), where we plotted F p (H). We 
observe that the kink in F p (indicated by the arrows) is still visible in the field range where the transition Hi should 
appear, and that it shifts to lower fields with increasing x, in agreement with the T dependence. In other words, the 
analogy between x and T is satisfied for H* , thus confirming that the kink arises from nonlocality. 

Another interesting fact is that the Lu-0 sample has a larger F p than the Y-0 sample at the same T, even though it 
has a lower density of magnetic impurities. This indicates that the magnetic moments in the Y-0 are not the relevant 
pins for the flux lines when H || c. 

B. Out-of-plane anisotropy 

We now turn to the pinning properties for H_Lc. Figure |^(a) shows F p at T = 3K for the Y-0 sample for H || [100], 
as a function of the reduced field h = H/H C 2- The H || [001] data, already shown in fig. |](a), is included for 
comparison. It is evident that the behavior for H || [100] is very different from that observed for H || [001]. First we 
note that a broad maximum develops at intermediate fields h max ~ 0.15. Second, F p is much larger than for H || c 
over most of the field range. The maximum of this out-of-plane anisotropy, T — F p [100]/ F p [001], is ~ 20 as shown in 
Figure ||(b). At higher T the overall behavior and the T values are similar. Finally, for H || [110] (not shown), we 
found that F p is slightly but systematically smaller (~ 12%) than for H || [100]. 



5 




FIG. 3. (a) Pinning force density for the samples Y-0 at T = 3K and (c) Lu-1.5 at T = 5K for H || c and H_Lc. (b) and (d) 
show the out of plane anisotropy T corresponding to (a) and (c) respectively. 

Hereafter we will focus on the origin of the difference in F p between the c-axis and the afe-plane. As we pointe d o u t 
in ref.cHl, the large out-of-plane F p anisotropy sharply contrasts with the very small (< 10%) mass aiiisotropycil'c3E3. 
Therefore, explanations based on the anisotropic scaling frequently used in high T c supexcpnductorsEj are ruled out. 

As we mentioned above, the Y-0 crystal contains a small amount of magnetic impuritiesEj, which align preferentially 
along the ab-plane and thus are a potential source of the F p anisotropy. However, we had previously arguedo that 
this was unlikely. Indeed, since the alignment of these localized moments increases with H, pinning should become 
more directional as the field increases and thus T should increase monotonically with H . This is not in agreement 
with the data in Figj^(b), where T first grows with h, maximizes at h ~ h max and decreases again. In particular, F p 
is almost isotropic at the peak effect (sec dotted line which corresponds to T = 1). Thus, although we were not able 
to totally rule out magnetic pinning on crystal Y-0, the data suggested that this was not the case. 

Conclusive evidence that the localized magnetic moments are not responsible for the large T comes from the 
persistence of this effect in the samples Lu-0, Lu-1.5 and Lu-3, which have a magnetic impurity content two order of 
magnitude lower than the Y-0 crystal. In Figure ||(c) we show F p (h) for the Lu-1.5 for H || c and H_Lc, at T = 5K, 
and in Figure. ||(d) we plotted the corresponding T(h). We observe that the out-of-plane anisotropy is even larger 
than in the Y-0 sample. In contrast to the behavior observed in the Y-0, T in the Lu-1.5 increases monotonically with 
h up to near the peak effect region, where it suddenly drops approaching to the isotropic limit. Measurements on the 
Lu-3 crystal at several T show a similar behavior and exhibit a T larger than in the Y-0. 

We had also previously arguedo that T in the Y-0 sample seems too large to be ascribed to nonlocality, which 
should appear as a perturbatively small effect. This conclusion is unambiguously confirmed by the presence of even 
larger T values in the doped samples Lu-1.5 and Lu-3, where nonlocal effects are strongly suppressed. 

The existence of significant surface barriers for H_Lc, have also been ruled out by performing minor hysteresis loops 
with H || ab at several T and H. Examples for the Y-0 and Lu-1.5 samples are shown in Figure [|. If hysteresis were 
due to surface barriers no flux changes would occur in the bulk while H is changing from one branch of the main loop 
to the other pije, hence the data of the minor loop connecting the lower and upper branches would be Meissner-like 
straight linesEH. In contrast, in the case of bulk pinning, the lines connecting both branches are curved (parabolic in 
the simplest Bean model for an infinite slab) just as we observe in the insets of Figure |(a) and (b). 
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H (kOe) 



FIG. 4. main panel: hysterisis loop for H_Lc at T = 5K for (a) Y-0 sample and (b) Lu-1.5 sample. The insets show a blow 
up of the minor hysterisis loops corresponding to the full loop showed in the main panel. The dotted lines show the behavior 
expected for surface barriers. 

Moreover, the F p calculated from these minor loops assuming only bulk pinning are in good agreement with those 
obtained from the main loops (Fig. Indeed, according to the isotropic Bean model the critical state profile can be 
completely inverted if we change the external field in AH = 2(4ir AM) / g , where g = 1 — t/3L, and t and L arc the 
shortest and largest sample dimensions perpendicular to the field direction, respectively. From the inset of figure ||(a) 
we find that the width of the hysterisis loop is AttAM « 150G. For this sample g ~ 0.9, thus it is expected that the 
connecting line between the two branches of the loop, reaches the upper branch, when the field is decreased about 
AH « 330G. Inspection of figure |](a) shows that AH w 500G, in reasonably good agreement given the simplicity of 
the model. The same analysis on the data of the inset of figure [I|(b), leds to a similar conclusion. 

IV. CONCLUSIONS 

Non-magnetic borocarbides provide a very attractive ground to study the limit of very low pinning forces. The 
high quality of the crystals have the double benefit of producing large mean free paths, resulting in measurable non- 
local effects, and a very dilute distribution of vortex pinning centers. The combination of both factors allows the 
observation of the influence of nonlocality on vortex pinning. This is a rather unexplored field, and further studies, 
both experimental and theoretical, arc needed. On the other hand, the origin of the large out-of-plane anisotropy 
in F p is unclear and deserves further investigation. We have conclusively ruled out explanations based on the mass 
anisotropy, pinning by magnetic impurities, non-local effects and surface barriers. A simple explanation for it could 
be the presence of some still unidentified anisotropic pinning centers, such as planar defects. 
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